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We consider the influence of quenched disorder on the relaxational critical dynamics of a system 
characterized by a non-conserved order parameter coupled to the diffusive dynamics of a conserved 
scalar density (model C). Disorder leads to model A critical dynamics in the asymptotics, however it 
is the effective critical behavior which is often observed in experiments and in computer simulations 
and this is described by the full set of dynamical equations of diluted model C. Indeed different 
scenarios of effective critical behavior are predicted. 
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The critical behavior of pure systems might be changed 
by introducing imperfections like dilution, defects, etc. 
into a critical system. If such a change can be expected 
is answered by the Harris criterion P| stating that a new 
diluted critical behavior appears if the specific heat of 
the pure system is diverging. The diluted critical be- 
havior then has a nondiverging specific heat. Since the 
borderline value Uc between a diverging and nondiverging 
specific heat at space dimensions d — 3 lies between or- 
der parameter (OP) dimensions n ~ 1 (Ising model) and 
n = 2 (XY model) only the Ising case belongs to a new 
universality class. In consequence this result led to the 
conclusion that for the critical dynamics the couplingof 
conserved quantities to the OP is of no relevance 0,131 ■ 
The argument was the following: For the critical dynam- 
ics of a relaxational model it was shown ^, IS that the 
coupling to a conserved density (e.g. the energy density) 
is relevant if the specific heat diverges. Due to dilution 
this is never the case and therefore the coupling is of no 
relevance. Therefore most of the papers considered only 
the relaxational dynamics of Ising systems ^, lOj 

However this argumentation is based on the asymp- 
totic properties of the diluted model. Experimental 
data and computer simulations made clear that in most 
cases one observes non-asymptotic critical behavior, 
described often by dilution dependent effective exponents 
(see e.g. JiJJJ). In such a case the Harris criterion does 
not hold and therefore one has to consider in the dynam- 
ics the coupling to the conserved density and its effects 
on the effective critical behavior. In addition one is not 
restricted to the Ising case since already in statics the 
effective critical behavior for n > 1 is different from the 



pure case |l2l |. 

There are two relevant parameters of model C: (i) the 
static coupling 7 of the OP to the conserved density and 
(ii) a dynamic parameter, the time scale ratio w = T/X 
where P is the relaxation rate of the OP and A is the 
diffusion rate of the conserved density. From the renor- 
malization group (RG) treatment of model C one knows 
that the one loop order does not give reliable results due 
to the stability of a fixed point with the time scale ratio 
w — 00. In two loop (and higher) order it turns out that 
this fixed point is unstable and model C is characterized 
by strong and weak scaling regions for the dynamics at 
d — 3 5, 6J. Moreover it was shown that non-asymptotic 
effects are already present in model C Q- In the follow- 
ing we will consider how these aspects are influenced by 
disorder. 

Model C 0, [3 describes the relaxational dynamics of 
a system characterized by an n-component nonconserved 
OP (po{x,t) coupled to the diffusive dynamics of a con- 
served scalar density mo{x, t). The structure of the equa- 
tions of motions is not changed by the presence of disor- 
der. They read: 
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where or ° denote unrenormalized quantities. The 
stochastic forces in ^ satisfy the Einstein relations: 

<e^^{x,t)0^^{x',t')> = 2tS{x~x')S{t-t')S^j, (2) 

<e,n{x,t)e,n{x',t')> = -2xy^s{x~x')S{t-t')(,3) 

Equilibrium is described by the static functional Ti. of 



FIG. 1: Two loop flow of static parameters for n = 1. The 
stable fixed point (full square) has a finite width A of the 
impurity distribution and leads to asymptotic exponents dif- 
ferent from the pure case (unstable fixed points are indicated 
by full circles). Depending on the initial values of the renor- 
malized couplings the effective exponents are quite different. 



the disordered magnetic system 
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where V{x) is an impurity potential which introduces 
disorder to the system, d is the spatial dimension. It 
contains a coupling 7 to the secondary density which 
can be integrated out. Thus static critical properties 
described by the functional Q are equivalent to those 

of the functional H = / (i'^x|if|(^oP + VixWo]"^ + 
\ Yll=i{'^^i,oY + ^l¥'o|'*|- The parameters f and u are 

related to f, u, am, 7 and hm hy r = f + jh,n/am and 
u = u — 37^/am. f is proportional to the temperature 
distance from the mean field critical temperature, u is 
positive. The properties of the random potential V{x) 
are governed by a Gaussian distribution with width A 
(<< V{x)V{x') »= 4:AS{x - x'), the double angular 
brackets means averaging over disorder). If 7 = equa- 
tions ^ describes dynamical properties of a purely re- 
laxational model (model A) in the presence of disorder 



We treat the critical dynamics of the disordered models 
within the field theoretical RG method 13], where the 
appropriate Lagrangians of the models are studied. The 
average over the random potential generates new terms 
in the Lagrangians with coupling A, which correspond to 
the static coupling terms in H generated by the disorder. 
The renormalization of the Lagrangian leads to the RG 
functions, describing the critical dynamics of our models. 
We use the minimal subtraction scheme with dimensional 
regularization to calculate these functions. 



FIG. 2: Two loop flow of static parameters for n = 3. The 
stable fixed point (full square) has a zero width A of the im- 
purity distribution and leads to asymptotic exponents equal 
to the pure case (unstable fixed points are indicated by full 
circles). However depending on the initial values of the renor- 
malized couplings the effective exponents are quite different 
and also depend on the inital value of the width A. 



For renormalization of the OP ipQ, fourth-order cou- 
plings u, A and correlation functions with 0^ insertion 

-« 1/2 

we introduce renormalization Z-factors as (po = Zip ip, 
u = fi'Z-^ZuuA^\ A = fi'Z-^ZAAA^\ and |</?oP = 
Z(p2|(/?p (/i is the scale, e — i — d and Ad is a geometric 
factor). Within dynamics renormalization factor for the 
OP kinetic coefficient F = Z^T is introduced. The Z- 
factors introduced so far are enough to renormalize the 
diluted model A. For model C one needs to introduce 
additional renormalization factors. The secondary den- 
sity mo and coupling parameter 7 are renormalized by 
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with 



Z-2(u, A, 7) = 1 + 7^A^2 (u, A). The kinetic coefficient 
A renormalizes as a„iA = Z^^^X. 

Defining the ^-functions as din Z^^/c?ln/i, where Z 
represents any renormalization factor, one obtains the 
flow equations for the renormalized static and dynamic 
parameters. The flow equations for u and A decouple 
from the remaining parameters and are equal to expres- 
sions obtained for any n in the diluted Ginsburg-Landau- 
Wilson (GLW) model JJJ. For the additional static pa- 
rameter 7 appearing in diluted model C we have 
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The flow parameter / is related to the reduced tempera- 
ture and is consequently a measure for the distance to the 
critical temperature. The function {u, A) is known 
from statics in the diluted model The function 

Bp2{u, A), which is defined by the additive renormaliza- 
tion of the specific heat within the GLW-model, in the 
diluted case reads Bp2{u, A) = n/2 + 0{u'^, A^, uA). 
The flow equation for the time scale ratio (we introduce 
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TABLE I: Static and dynamic fixed points for n = 1. Fixed 
points witfi p* = 1 (ui* = oo) axe always unstable and not 
shown. M is the stable fixed point of the diluted model C. 
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The corresponding explicit two loop expression for the C- 
function of model C, C,f^\u,^, p) = (^\u,^, p) +nj'^ /2, 
is given in P| (see Cw in Eq. (50) there). In these terms 
also the pure model A terms are included, which taken 
together with the last three terms in the first line of Eq. 
O recover the (- function for the diluted model A 0, U ■ 
The zeros of the right hand sides of Eqs.©, and the 
corresponding equations for u and A give the possible 
fixed points. 

It turns out that for all values of n the stable fixed 
point value 7* = 0, and the well-known asymptotic static 
results are reproduced. Thus the flow in the space of the 
static couplings u, A and 7 for n < < 2 looks like the 
flow for n = 1 (see Fig. ^ whereas for n > ric it looks 
like the flow for n = 3 (see Fig. The fixed points 
for n — 1 are indicated in Tab. ^ Only the mixed fixed 
point (M in Fig^ is stable. However depending on the 
initial conditions a rich crossover behavior is observed. 



TABLE II: Stability exponents according to the fixed points 
in TabElfor n = 1. 
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FIG. 3: Effective dynamical exponent for n= \ for different 
initial conditions. The lines correspond to the static flows 
shown in Fig0 

The same is true for n = 3 although now the pure fixed 
point (P in Fig. [SJ is stable. 

The stable fixed points are found by calculating the 
stability matrix and its eigenvalues loi {i represents it, A, 
7 or p) . They govern the "velocity" of the flow near 
the fixed points. A small stability exponent at the stable 
fixed point indicates a slow approach of the asymptotic 
behavior. This is the case for all n and can be seen from 
Tab im for n — \. It is the slow approach in the 7 direction 
which characterizes the static flow within model C (see 
the extremely small values lo^ = 0.0018 for n = 1 [at the 
fixed point M in TabHIl and torn = 3 stiU w-^ = 0.1109). 
In addition to the small values of the static stability ex- 
ponents near the fixed points one has also small values 
ijOp coming from the dynamic parameter p (however only 
near unstable fixed points). Thus one expects depending 
on the initial values of static and dynamic parameters a 
complex behavior in the non-asymptotic region. The 
effective dynamical exponent is of special interest here, 
it is found by inserting the solutions of the flow equations 
into the expression 

Ze//(0 = 2 + Cr(M(0,A(0,7(0,P(0)- (8) 

One reaches the universal asymptotic value z when the 
flow comes very near the stable fixed point. One observes 
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FIG. 4: Effective dynamical exponent for n = 3 for different 
initial conditions. The lines correspond to the static flows 
shown in Fig0 
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FIG. 5: Different contributions to the effective dynamical 
critical exponent for n = 1. Tfie solid line represents the 
total Zeff, for the other lines see the text. 

that the stable fixed point value of the time ratio is always 
zero, independently of the specific heat exponent value of 
the pure model. Consequently, the results for model A 
are recovered in the asymptotics (either in the diluted 
universality class for n < or in the pure model A 
universality class for n > ric). This combines with the 
fact that even in the region where dilution changes the 
static critical behavior the stable fixed point value of 7* 
is always zero. In consequence the secondary density is 
for all n asymptotically decoupled and has an asymptotic 
dynamical exponent = 2. 

The non-asymptotic behavior is however quite differ- 
ent as can be seen from Figs. Inland 01 Due to the static 
non-asymptotic behavior also the dynamics is dominated 
by different non-asymptotic effects. This can be seen 
by comparing different 2e//(0 different initial condi- 
tions. The curves a, b, and c in Fig. ^ ^-^d curves 
in Fig. 121 reach large values of the coupling 7 and/or A 
and this leads to the typical maximum in the effective 
exponents independent of the initial value of p (for the 
statics see e.g. ^3)- However an additional fixed point 



P' is present at n = 1. This leads for curve d in Fig. 
almost to a plateau of z^f / at its value for the unstable 
fixed point P'. This plateau is more pronounced when 
the flow comes nearer to P' where it stays longer because 
of the small transient exponent ujp. For curve c both ef- 
fects (the maximum and the effect of fixed point P') are 
combined leading to the minimum in Zefj. Consider now 
the contributions to the effective dynamical critical expo- 
nent Zef f of different origin: (i) from the terms already 
present in model A (dashed curve in Fig. IH)), (ii) from 
terms present in pure model C only (short dashed curve) 
and (iii) and from terms present in the diluted model 
C only (dashed-short dashed curve). The above contri- 
butions may add up to almost the asymptotic value 
of the exponent although the parameters are far away 
form their asymptotic values. This is an important point 
since the appearance of an asymptotic value in one phys- 
ical quantity does not mean that other quantities have 
also reached the asymptotics. This is due to the different 
dependence of physical quantities on the model param- 
eters. Another special feature of the diluted model C 
is that already in one loop order one observes qualita- 
tively the same behavior as in two loop order of course 
with changed values for the exponents and the borderline 
value Tie, which in one loop is at Uc — 4. 

In concluding we remark that contrary to the general 
belief the coupling of a conserved density to the order pa- 
rameter is relevant for the calculation of the dynamical 
critical behavior of diluted systems since this coupling 
is important to describe non-asymptotic effects. These 
effects have been seen in the experiments on physical 
systems fl3. (T^j l as well as in Monte Carlo simulations 
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